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Abstract 

Numerical solutions of the standard electrokinetic model provide a basis for inter- 
preting a variety of electrokinetic phenomena involving 'bare' colloids. However, the 
model rests on the classical notion of a shear or slipping plane, whose location is 
unknown when surfaces are coated with permeable polymer. Consequently, an elec- 
trokinetic model for 'soft', 'hairy' or 'fuzzy' colloids has been developed, but until 
recently solutions were available only for several restricted cases, most notably for 
particles with thin, uniform layers, and without polarization and relaxation. Here 
we present numerically exact solutions of the full model for a variety of soft col- 
loids, including PEG-coated liposomes, PEO-coated latices, human erythrocytes, 
and polyelectrolyte micelles. Particular attention is given to linking the thickness, 
density and permeability of the coatings, which are key parameters in the model, 
to "physical" quantities, such as the polymer molecular weight, adsorbed amount, 
and hydrodynamic layer thickness. This paper also identifies limits — on the ionic 
strength, particle size, layer thickness and permeability — beyond which earlier the- 
ories breakdown. In short, polarization and relaxation are as influential on the mo- 
bility of soft colloids as they are for 'bare' particles. 

Key words: electrophoretic mobility, electrokinetic phenomena, soft colloids, 
polymer-coated colloids, sterically stabilized colloids, stealth liposomes, human 
erythrocytes, charged micelles 
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1 Introduction 



Polymer deposited on the surface of 'bare' colloidal particles can be used to 
regulate the stability and rheology of dispersions; the influence of thin polymer 
coatings on the drag coefficient and intrinsic viscosity, for example, are avail- 
able from exact solutions of the Stokes equations (jlj). In addition, adjusting 
the electrostatic, steric, and dispersive contributions to the inter-particle po- 
tential affords control over the kinetics and thermodynamics of self-assembly 
processes. In turn, these bear directly on the microstructure of materials with 
novel properties. Nevertheless, the polymer coating density, thickness and 
charge affect electrokinetic transport processes near the supporting surface, 
making it difficult to discern competing influences that affect light scatter- 
ing and electrophoresis measurements, for example. Clearly, the application of 
polymer coatings for tailoring the properties of dispersions will benefit from a 
reliable electrokinetic theory. 

The purpose of this paper is to describe applications of a recently developed 
model for predicting the influence of polymer coatings on a variety of elec- 
trokinetic particle characterization measurements. The focus of this paper is 
on the electrophoretic mobility, but it should be kept in mind that the model 
permits solutions under dynamic conditions, as encountered in dielectric re- 
laxation spectroscopy and electro-sonic-amplitude devices. 

While the electrokinetics of bare particles are (relatively) well understood, the 
situation for colloids coated with a permeable polymer is less established. With 
bare particles, ion transport and hydrodynamics follow classical models, where, 
for example, hydrodynamic shear begins at the particle surface (shear surface). 
Efforts to model transport processes inside a thin porous layer have adopted a 
variety of approaches (0;0)- More realistic (and physically appealing) models 
treat polymer layers as continuous porous media with a permeability that 
reflects the density and hydrodynamic size of the polymer segments. 

Charged layers, which are important because of their biological context, present 
a challenging set of problems. A variety of approximations have been used to 
construct tractable models, but, until recently, none had the generality of 
O'Brien and White's numerical solution of the standard model for bare par- 
ticles (J4|). For example, Brooks and coworkers developed models for human 
erythrocytes (red blood cells), with 'flat', uniformly charged layers to mimic 
the extracellular dvcocaylx layer fii); Hermans and Fujita considered porous 
charged spheres and , more recently, Ohshima and coworkers general- 

ized these models, presenting approximate solutions (for the electrophoretic 
mobility) that neglect, in particular, polarization and relaxation (fil; liof). 

Without exact solutions, the accuracy of approximations is unknown, so empir- 
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ical parameters derived from experiments using approximate theories should 
be accepted with caution. Saville demonstrated that even 'large' particles with 
thin polymer and diffuse double layers are susceptible to polarization and re- 
laxation Inferences drawn from his semi-analytical theory have been 
borne by numerically exact solutions of the full model fl^). and recent atten- 
tion has turned to details of the polymer segment density distribution ffl3h . 



This paper reviews the full electrokinetic model for colloids with soft neutral 
and charged layers. We draw upon numerically exact solutions of the gov- 
erning equations, as provided by Hill, Saville and Russel's model (H2). Their 
methodology removes all restrictions imposed by earlier approximate theories, 
permitting studies that examine the influence of particle size and polymer- 
layer structure in considerable detail. Accordingly, this paper establishes the 
parameter spaces of earlier approximate theories as limiting cases of the full 
model. 



We begin with a presentation of the full electrokinetic model (§ 2), linking the 
solution of the equations to the mobility and effective layer thickness (§ 2.7), 
and briefly discussing key parameters that characterize polymer layers (§ 2.8). 
Next, the results are presented, beginning with particles that have neutral 
polymer coatings (§ 3). First, we highlight the significant role of the hydro- 
dynamic layer thickness, distinguishing it from the actual thickness (§3.2). 
Detailed comparisons are made with Ohshima's theories for thin coatings at 
high ionic strength, and, in particular, establishing the significance of polariza- 
tion and relaxation (§ 3.3). Next, a methodology for interpreting experimental 
data is presented (§ 3.4), and the full model is then applied to interpret the 
mobilities of stealth liposomes (§ 3.5) and PEO-coated latices (§ 3.6). Turning 
to particles with charged layers (§ 4), we compare the full model with ap- 
proximate theories for thin layers (§ 4.1). 'Exact' calculations are presented 
for human erythrocytes (§ 4.2) and 'small', highly charged polyelectrolyte mi- 
celles (§ 4.3). We conclude with summary (§ 5). 



2 Theory 



The electrokinetic model and the computational methodology adopted in this 
work have been described in detail by Hill, Saville and Russeljl2|). The model 
augments the standard electrokinetic model for bare colloids (|4j) with arbitrary 
(radial) distributions of Stokes resistance centers and charge. Electroosmotic 
flow within the polymer layer is hindered by hydrodynamic drag on the seg- 
ments. In this work, elect romigr at ion and molecular diffusion are unaltered, 
which is reasonable when (i) the hydrodynamic volume fraction of the ions is 
low and (ii) the ions are small compared to the polymer intersties, i.e., the 
polymer volume fraction is low. In principle, however, the methodology per- 
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mits the introduction of spatially dependent ion mobilities, solvent viscosity 
and dielectric constant. 

2.1 The electrokinetic transport equations 

The transport equations and boundary conditions are presented here in di- 
mensional form. As usual, they comprise the non-linear Poisson-Boltzmann 
equation 

N 

e e s V 2 ^ = - J2( n J ~ n j) z j e ' (!) 
j'=i 

where e G and e s are the permittivity of a vacuum and dielectric constant of 
the electrolyte (solvent); rij and are the concentrations of the jth mobile 
and fixed charges, with Zj the valences 1 ; e is the elementary charge and ip the 
electrostatic potential. 

Transport of the mobile ions is governed by 

Qnriajiu- vj) - ZjeVtp - kTVlnnj =0 (j — 1, N), (2) 

where a,j are Stokes radii, obtained from limiting conductances or diffusivities; 
rj is the electrolyte viscosity, and u and Vj are the fluid and ion velocities; kT 
is the thermal energy. 

Ion conservation demands 

dn j /dt = -V-(n j v j ) (j = l,...,A0, (3) 

where t denotes time, with the ion fluxes j ■ = rijVj obtained from Eqn. (2). 

Similarly, momentum and mass conservation require 

N 

p s du/dt = r]V 2 u - Vp + r]/£ 2 (V - u) - ^ njZjeVtp (4) 

i=i 

and 

V-u = 0, (5) 

where p s and u are the electrolyte (solvent) density and velocity, and p is 
the pressure. Note that r]/£ 2 (V — u) represents the hydrodynamic drag force 
exerted by the polymer on the interstitial fluid, with £ 2 the (radially varying) 
permeability of the polymer, and V is the particle velocity. Clearly, the poly- 
mer and particle are assumed to move as a rigid composite, and the particle 
Reynolds number Re = Vap s /i] is assumed small. 

1 The valence of the fixed charge is set opposite to that of its respective (mobile) 
counterion in Eqn. (1) 
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As usual, the double-layer thickness 



k- 1 = [kTe s e /(2Ie 2 )] 1 / 2 (6) 
emerges from Eqns. (1) and (2) where 

N 

/= a/2) E^r (7) 

is the bulk ionic strength, with n°° the bulk ion concentrations. Note that ion 
diffusion coefficients are 

Dj = kT / (Qnrjaj), (8) 

and the permeability (square of Brinkman screening length) may be expressed 
as 

i 2 = \/{n^a s F s ) = 2a 2 /(90 s F s ), (9) 

where n(r) is the concentration of Stokes resistance centers (segments), with 
a s and F s (<j) s ) the Stokes radius and drag coefficient of the segments. For the 
rigid, 'random' configurations of spherical segments assumed in this work, 
the segment drag coefficient varies with the hydrodynamic volume fraction 
4> s = n(A/3)nal according to 

= f + 3(0,/2) 1 / 2 + (135/64)0, In S + f 6.4560, 

f + 0.6810, - 8.4802 + 8.1603 {<Ps ' 1 1 ' 



The numerator in Eqn. (10) derives from theory for low to moderate volume 
fractions (flit), with the denominator providing a smooth transition through 
data from multipole simulations (0, < 0.4) to the well-known Carman corre- 
lation (0.4 < 0, < 0.64). Note that other drag-coefficient correlations exist for 
porous media with ordered and other random microstructures (see 2 12, for 



fibrous microstructures). At present, however, such details are of secondary 
importance to the distribution of segments n(r) and their size a s . The seg- 
ment and fixed-charge density distributions are discussed in section 2.8, and 
throughout sections 3 and 4. 



2.2 Inner (particle surface) boundary conditions 



The model permits either the equilibrium surface potential ( or surface charge 
density a c to be specified. Because the surface (r = a) is assumed impermeable, 
the surface charge is constant under the influence of an external electric field 
or particle motion. This is achieved with no-flux boundary conditions for each 
(mobile) ion species. As usual, the no-slip boundary condition applies to the 
electrolyte. It follows that (inner) boundary conditions are either 

if) = C at r = a (11) 
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or 

e s e Vtp\ 0Ut ■ h - e p e Vtp\ in ■ n = -a c at r = a, (12) 

with 

rijVj ■ h = at r = a (j = 1, N) (13) 

and 

u = V at r = a. (14) 

Note that n = e r is an outward unit normal and r = re r is the (radial) 
distance from the center of the particle. The dielectric constant of the bare 
particle e p does not affect the electrophoretic mobility to linear order in the 
applied electric field (0), but it is relevant under dynamic (high-frequency) 
forcing, such as in dielectric spectroscopy. 



2.3 Outer (far-field) boundary conditions 

Particle interactions are neglected, so an otherwise stationary, infinite domain 
leads to far-field boundary conditions 

ip — > — E ■ r as r — ► oo, (15) 
rij — > n°° as r — > oo (j = l,...,N), (16) 

and 

u — >• as r — >• oo, (17) 
where -E is a uniform applied electric field. 

2.4 Solution of the equations 



Coupling the transport equations and boundary conditions to the particle 
equation of motion specifies the electrokinetic model 2 . As usual, the solution 
requires non-linear perturbations to the equilibrium state to be neglected. 
Perturbations (proportional to the applied electric field) are introduced via 

if) = V>° -E-r + ip', (18) 



n,=n° + n; (j = 1, N), (19) 

and 

V = P° + p', (20) 



2 A computer package (MPEK) implementing Hill, Saville and Russel's methodol- 
ogy is currently available from the corresponding author. 
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where the first terms on the right-hand sides denote the equilibrium values, 
and the primed quantities denote the perturbations. 



2.5 Equilibrium state 



When E = V = 0, equilibrium is specified according to 

N 

e e s VV° = - EK - ^)% e > (21) 
j'=i 

= _V • [-Dj VnJ - Z] eD 3 /{kT)ny^} (j = 1, TV), (22) 

and 

= - Vp° - ]T n° jZj eV^°, (23) 
i=i 

with either 

^° = C (24) 

or 

e s e Vip°\ out • e r - e p e VV>°L • e r = -cr c at r = a, (25) 



[-DjVnJ - z j eD j /(kfT)n^ip°] ■ e r = at r = a (j = 1, TV), (26) 

and 

^° -> as r -> oo, (27) 

n° -> raf as r -> oo (j = 1, N). (28) 

2.6 Linearized perturbed state 
Under periodic forcing 

E = Eexp (-iut)e z , (29) 
with frequency u/(2tt), axisymmetric perturbations take the forms 

rj}' = j>(r)E ■ e r , (30) 

n'j = n 3 {r)E ■ e r (j = l,...,N) (31) 

and 

u = -2(h r /r)(E ■ e r )e r - (h rr + h r /r)(E ■ e e )e e , (32) 
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where 

u = V x V x h(r)E. (33) 

The perturbations satisfy 

N 

e e s VV = -5>;-z,e, (34) 
i=i 

dn , j /dt=-V-j' j (j = l,...,iV), (35) 

where 

f = -DjVn'- - zjeDj/ikTy-Vip - ZjeDj/ikT^AVtp' - E) (36) 



N 

p s du/dt = T]V 2 u - Vp' - r]/£ 2 (u -V)-J2 n° jZj e(Vtp' - E) (37) 

AT 



and 
with 



and 



V • u = 0, (38) 

e s e (V^' - #)U* • e r - e p e (V^' - E)\ in ■ e r = at r = a, (39) 

j-;- e r = at r = a (j = 1, N), (40) 

u = y at r = a, (41) 

^' -> as r -> oo, (42) 

nj . -> as r -> oo (j = 1, N), (43) 

u — > as r — > oo. (44) 

,2. 7 Measurable quantities derived from the solution 



The linearized equations and boundary conditions, which are characterized 
by multiple and widely varying length scales, are solved using the methodol- 
ogy described in detail by Hill, Saville and Russel fsee Il2h . In practice, the 
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equations are solved with the particle fixed at r = 0, with (i) a prescribed 
far-field velocity —V and E = — the so-called (U) problem, or (ii) a pre- 
scribed electric field E and stationary far-field — the so-called (E) problem. 
These solutions are superposed to satisfy the particle equation of motion (see 
S EH), which leads to the electrophoretic mobility, 

V/E = C E /C U } (45) 

where 



C x = lim hf. (46) 

r— >oo 

Under oscillatory conditions, the dynamic mobility is 

V/E = C E /[C U - M 3 ( Pp - p s )/(3p s )], (47) 

where 

C x = lim hfr 2 (48) 

r— >oo 

and p p is the density of the bare particle. Note that the superscripts X e 
{[/, .E} denote the (U) and (E) problems, where C x are "asymptotic coeffi- 
cients" that characterize the far- field decay of h(r) in Eqn. (33). 

Following standard convention, the scaled electrophoretic mobility is defined 
as 

M = (V/E)3 V e/(2e s e kT), (49) 

so the well-known Hiickel and Smoluchowski mobilities are M = ( and M = 
(3/2)£, where ( is scaled with kT/e. 

The drag coefficient 

F = f/(6irr}aV) (50) 

is the force / required to translate the particle with velocity V in the absence 
of an applied electric field, scaled with the Stokes drag force on the bare 
uncharged sphere. Frictional drag on the polymer coating and electro-viscous 
effects arising from charge increase the effective particle size. This motivates 
our definition of an effective coating thickness 

L e = a(F-l), (51) 

which, in the absence of electro-viscous stresses, is termed the hydrodynamic 
coating thickness 

Lh = L e when a c = nj = 0. (52) 

In practice, the hydrodynamic coating thickness is realized when the surface 
potential is low, which is often the case at high ionic strength. 
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2.8 Polymer layer characterization 



The length and molecular weight of a polymer segment, / and M s , and the 
number of segments per chain N are defined so the chain contour length and 
molecular weight are l c = Nl and W = NM S . If the segments are monomer 
units, then I and N are denoted l m and N m , so l c = N m l m and W = N m M m . 
Similarly, if the segments are defined as statistical or Kuhn segments, then 
the mean-squared end-to-end distance of a free chain in a theta solvent is 
R = IN 1 / 2 = J[(M/M TO )(W0] 1/2 > g^ing 

l = R 2 M m /{Ml m ). (53) 

This is useful because R can be expressed in terms of the hydrodynamic radius 
and radius of gyration of 'free' polymer chains, and l m is usually known from X- 
ray crystallography or knowledge of bond lengths and monomer architecture. 

The mass of polymer per unit area of substrate is 

Wo p = NM s a p = N m M m a p , (54) 

where a p is the number of chains per unit of surface area. The number of 
segments per unit volume n is often scaled to give a volume fraction <fi = 
nl 3 . The 'actual' volumes of a solvent molecule and a polymer segment are 
v\ and t>2, so the physical volume fraction of polymer is m>2- Similarly, the 
hydrodynamic volume fraction of polymer segments is <ft s = n(4/3)iral, where 
a s is the hydrodynamic (Stokes) radius of a segment. 

The grafting density is often scaled to give a p l 2 , but a better qualitative 'mea- 
sure' of the grafting density is, perhaps, a p R 2 = a p Nl 2 = a p {W / M m )l m l. 
When a p Nl 2 ~ 1, the layer is considered semi-dilute < 1) and, there- 
fore, laterally homogeneous. When a p Nl 2 < 1, however, the layer is dilute 
and, therefore, laterally inhomogeneous. Moreover, when a p Nl 2 > 1, layers 
are homogeneous throughout and are referred to as brush-like. This 'picture' 
applies to di-block copolymers where one (shorter) block anchors the other 
(longer) block to the substrate. The segments of adsorbed homopolymers in- 
teract strongly with the surface, giving rise to a dense inner region with a 
sparse periphery. The structure of polyelectrolytes is more complicated, since 
the statistical segment length varies with the bulk ionic strength and the de- 
gree of counterion dissociation, the latter of which may depend on the ionic 
strength and pH flih. 
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3 The mobility of colloids with neutral coatings 

3.1 Connection of the electrophoretic mobility to the hydrodynamic layer 
thickness 

In contrast to approximate analytical solutions, numerically exact solutions 
of the full electrokinetic model permit arbitrary segment and charge density 
distributions. Here we compare the electrophoretic mobilities of particles with 
uniform and non-uniform (Gaussian-like) layers, each with the same adsorbed 
amount and the same hydrodynamic layer thickness. In turn, the hydrody- 
namic layer thickness emerges as a characteristic of prime importance. 

Gaussian-like distributions 

n(r) = ni exp [— (r — a) 2 / 5 2 } (r > a), (55) 

where n\ specifies the density at the bare surface (r = a) and 5 characterizes 
the layer thickness, mimic terminally anchored polymer in a theta solvent when 
the molecular weight and grafting density are moderately low. Substituting 
Eqn. (55) into the general relationship 

/•OD 

Na p a 2 = / n(r)r 2 dr, (56) 

J a 

where a p is the chain grafting density and iV is the number of segments per 
chain, gives 

Na p a 2 = (n 1 5 3 /A){ v ^[2{a/5) 2 + 1] + Aa/5}. (57) 
Similarly, for a uniform distribution with density n and thickness L, 

Na p a 2 = (n /3)a 3 [(l + L/af - 1]. (58) 

Let us consider how either n\ and 5 or uq and L are related to the elec- 
trophoretic mobility and hydrodynamic layer thickness Lh, both of which are 
routinely measured. Because Lh depends on the permeability £ 2 = l/(67in(r)a s ), 
the Stokes radius of the segments a s must be correctly specified; here, a s = 
0.175 A is a value for poly(ethylene glycol) (PEG) that provides a good 'fit' 
of theory (with a self-consistent mean-field description of the polymer density 
distribution) to measured electrophoretic mobilities (13). The physical char- 
acteristics of four representative coatings, each decorating the surface of an 
impermeable spherical colloid with surface charge density o c pa —0.78 /xCcm -2 
and radius a = 500 nm are listed in table 1. Each of the four sets of parameters 
(cases 1-4), which specify and iV (columns 2 and 3), have been used to char- 
acterize Gaussian-like and uniform (step-like) segment density distributions. 
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Table 1 



Parameters that characterize the layers whose segment density distributions and 
electrophoretic mobilities are shown in figures 1 and 2: a p l 2 = 0.072 with I = 
0.71 nm; a s = 0.175 A; a = 500 nm; a c = -0.780 ^Ccm' 2 . 









Gaussian-like 




Step-like 




case 


L h (nm) 


TV 


m (M) 


S (nm) 


n (M) 


L (nm) 


I (nm) 


1 


28 


46.2 


0.63 


18.8 


0.32 


31.8 


3.9 


2 


23 


18.5 


0.25 


18.8 


0.14 


29.3 


6.0 


3 


14 


18.2 


0.38 


12.5 


0.22 


18.7 


4.8 


4 


7.4 


9.10 


0.25 


9.4 


0.16 


12.8 


5.6 




scaled distance, (r-a)/a 



Fig. 1. The (scaled) segment density distributions <j){r) = n(r)l 3 as a function 
of the (scaled) distance from the bare surface (r — a) /a for the particles whose 
electrophoretic mobilities are shown in figure 2. The density and thickness of the 
uniform layers {dashed lines) have been adjusted to have (i) the same hydrodynamic 
layer thickness and (ii) the same number of polymer segments Na p as the 
corresponding non-uniform layer. See table 1 for the coating parameters, and figure 2 
for the electrophoretic mobilities. 

Note that the hydrodynamic layer thicknesses were established iteratively from 
solutions of problem (U) at high ionic strength. Clearly, the various molecular 
weights and hydrodynamic coating thicknesses lead to a variety of physical 
thicknesses (5 and L) and densities (no and ni). As shown in figure 1, the 
extent (maximum density) of each non-uniform layer is significantly greater 
(less) than that of its uniform counterpart. 

The electrophoretic mobilities are shown in figure 2 over a range of ionic 
strengths. Despite large variations in density and extent, constraining each 
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10 1 io 2 io 3 lu io' 10" 10" 

scaled reciprocal double-layer thickness, kw scaled reciprocal double-layer thickness, Ka 

(a) (b) 

Fig. 2. The (scaled) electrophoretic mobility M = 3r]eV/ (2e s e kTE) as a function 
of the (scaled) reciprocal double-layer thickness Ka (aqueous NaCl at T = 25° C with 
radius a = 500 nm). The surface charge density o c ~ —0.780 /^Ccm -2 at all ionic 
strengths. Solid lines are the full model with non-uniform segment density distribu- 
tions (solid lines in figure 1), dashed lines are the full model with uniform (step-like) 
segment density distributions (dashed lines in figure 1), and dash-dotted lines are 
Ohshima's theory Eqn. (11.4.27)). See table 1 for the coating parameters (cases 
1-4, with case indicating the bare particle), and figure 1 for the segment density 
distributions. 

uniform layer (solid lines) to have the same hydrodynamic thickness as its 
respective non- uniform layer (dashed lines) leads to similar mobilities. This is 
to be expected at low ionic strength, but the close correspondence at high ionic 
strength is surprising. When the diffuse double layer resides mostly within the 
layers (at high ionic strength), the mobility is expected to be more sensitive 
to details of the distribution. However, when the densities and thicknesses are 
constrained by Eqn. (56), the mobility reflects — to a good approximation — the 
hydrodynamic layer thickness. Clearly, if the hydrodynamic layer thickness is 
known, the electrophoretic mobility has a remarkably close connection to the 
underlying surface charge over a wide range of ionic strengths. 

The logarithmically scaled mobility axis (panel (b) in figure 2) facilitates a 
closer comparison of the numerically exact solutions with Ohshima's well- 
known analytical approximation for thin uniform layers (0 Eqn. (11.4.27)) 
(dash-dotted lines). Good agreement is seen when L/a < 0.02 and Ka > 100, 
but, as expected, the approximate theory breaks down at low ionic strength. 
As shown below, Ohshima's more complicated formula (19, Eqn. (11.4.24)) 
can be advantageous with higher electrostatic potentials. Here, however, the 
surface charge density is low enough for the breakdown to be attributed to 
the finite particle size, i.e., to polarization and relaxation. 
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Fig. 3. The (scaled) electrophoretic mobility M = 3r]eV/(2e s e kTE) as a function 
of the (scaled) reciprocal double-layer thickness kcl (aqueous NaCl at T = 25° C 
with radius a = 1750 nm). The surface charge density o~ c ~ —1.95 //Ccm -2 at 
all ionic strengths. Solid lines are the full model with uniform segment den sity 
distributions, and the dashed and dash-dotted lines are Ohshima's theory for low (|19l . 
Eqn. (11.4.27)) and arbitrary ijliL Eqn. (11.4.24)) electrostatic potentials. See table 2 
for the coating parameters. 

3.2 Thin coatings and passage to the high-ionic- strength limit 



Surprisingly, the full model deviates from Ohshima's flat-plate theories at 
high ionic strength, even with very thin layers. For example, figure 3 shows 
the electrophoretic mobilities of particles whose (uniform) coatings are all very 
thin compared to the particle radius (a = 1750 nm). The coating parameters 
and various dimensionless ratios are listed in table 2. Note that all the layers 
(cases 1-6) have the same grafting density and polymer molecular weight, 
so the increase in thickness L is accompanied by a decrease in density no 
and, hence, an increase in permeability £ 2 . Nevertheless, Brinkman screening 
ensures that the hydrodynamic layer thickness increases with L. 

Despite a substantial decrease in L/a — relative to the cases presented in fig- 
ure 2 — there is a significant difference between the approximate and 'exact' 
results when L/a > 0.01. Remarkably, at high ionic strength, the mobility 
decreases slowly with decreasing layer thickness, whereas the flat-plate theo- 
ries (dashed and dash-dotted lines) infer a monotonic decline, yielding much 
lower mobilities. We attribute this unexpected behavior to the finite particle 
size and a transition to the regime where kI 3> 1 with l/L <C 1. Note that 
the pressure gradient must be identically zero with a perfectly flat interface. 
However, when kL ^> 1, fluid in the electrically neutral region of the polymer 
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Table 2 

Parameters that characterize the uniform layers on particles whose electrophoretic 
mobilities are shown in fig ure 3: a p l 2 = 0.072 with I = 0.71 nm; N = 100; a, = 
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case 


Lh (nm) 

ft, \ ) 


n (M) 


L (nm) 


I (nm) 


£/L 

1 


LI a 

1 


1 


1.10 


14.0 


1.69 


0.591 


0.35 


6.3 x 10~ 4 


2 


4.37 


4.34 


5.45 


1.07 


0.20 


2.7 x 10~ 3 


3 


8.75 


2.31 


10.2 


1.47 


0.14 


0.005 


4 


17.5 


1.20 


19.6 


2.04 


0.10 


0.01 


5 


35.0 


0.613 


37.9 


2.86 


0.075 


0.02 


6 


70.0 


0.307 


74.0 


4.04 


0.055 


0.04 



(n^ 1 < r — a < L) must be driven by viscous stresses when k£ > 1. With a fi- 
nite radius of curvature, radial and tangential pressure gradients then develop 
to move fluid through the polymer, into and out of the underlying diffuse dou- 
ble layer. As demonstrated by Hill, Saville and Russel (flii Fig. 6), decreasing 
the permeability increases the tendency of fluid to enter and exit the poly- 
mer layer radially, and this, presumably, decreases the radial and tangential 
pressure fluctuations. 

As shown in figure 3, the breakdown of the flat-plate approximation at high 
ionic strength can be significant for neutral layers with moderate thickness and 
permeability. For example, the mobility with L w 20 nm and I ~ 2 nm (case 
5) deviates from Ohshima's theory when the ionic strength exceeds 18 mM 
(na > 770). Note that there exists an intermediate range of ionic strengths, 
with kI sufficiently small and na sufficiently large, where Ohshima's theory 
provides an excellent approximation. As expected, the higher surface charge 
here — relative to the cases presented in figure 2 — demonstrates that Ohshima's 
formula {HJ Eqn. (11.4.24)) (dash- dotted lines) for arbitrary potentials can in- 
deed provide a significantly better approximation than his simpler formula (flil 
Eqn. (11.4.27)) for low potentials (dashed lines). 

3.3 Polarization and relaxation, and the influence of particle size 

With polymer layers, clearly identifying the regions of the parameter space 
where polarization and relaxation are significant is not straightforward. The 
calculations presented in figures 4 and 5 systematically explore the influence of 
particle size (a ~ 50-3500 nm) and layer thickness. Similarly to figure 3 with 
a = 1750 nm, all the coatings have the same grafting density (a p l 2 = 0.072) 
and polymer molecular weight (N = 100), so the permeabilities vary with the 
(specified) hydrodynamic thickness and particle radius. The numerically exact 



16 



results (solid lines), which are compared to Ohshima's theory for arbitrary 
surface potentials ()19l Eqn. (11.4.24)) (dash- dotted lines), demonstrate how 
particle size influences the mobility. Because the surface charge densities are 
the same, varying the surface curvature is accompanied by a change in the 
surface potential at fixed ionic strength (k ~ y/l). Increasing the coating 
thickness decreases the mobility, mimicking the effect of decreasing the charge. 
As expected for bare particles, decreasing the radius increases polarization, 
which, in turn, decreases the mobility and shifts the maxima to smaller values 
of Ka. 



Neutral coatings attenuate polarization, however, as demonstrated by the in- 
creasing accuracy of Ohshima's theory with increasing layer thickness when 
Ka ~ 300 and a = 3500 nm (last panel in figure 5). The close correspon- 
dence for the thickest coating (L/a ~ 0.02) establishes the validity of the 
flat-plate approximation and the negligible influence of polarization under 
these conditions. It follows that decreasing the layer thickness should improve 
the flat-plate approximation. However, the agreement clearly diminishes, and 
this must be attributed to polarization and relaxation. Note that decreasing 
the layer thickness (with na = 300 and a = 3500 nm) yields kL = 1 when 
L pa 12 nm. Therefore, the four thinnest layers have kL < 1 with Ka > 1, this 
being the regime in which Saville ([lit ) developed a semi-analytical theory to 
highlight the role of polarization and relaxation. Comparing the exact results 
with Ohshima's theory (last panel in figure 5) confirms that polarization and 
relaxation can be significant when na is large. Clearly, when most of the double 
layer resides beyond the polymer, i.e., when kL <C 1, the ratio of polarization 
(by convection and electromigration) to relaxation (by molecular diffusion) is 
expected to be maximal. 



3.4 Interpreting experiments 



As shown in section 3.1, the full electrokinetic model provides satisfactory 
predictions of the electrophoretic mobility with uniform layers when the coat- 
ing density and thickness reflect the correct adsorbed amount a p N (according 
to Eqn. (58)) and the 'correct' hydrodynamic thickness Lh- This approach 
requires the Stokes radius of the polymer segments a s to be specified. Note 
that, in the Brinkman model, the designation of monomer or statistical units 
to a polymer segment is arbitrary, as long as (i) the number of segments is un- 
ambiguously related to the molecular weight, e.g., W = NM S = N m M m with 
l c — Nl — N m l m , and (ii) the Stokes radius of the segments is appropriate. 

In general, determining the hydrodynamic thickness of a layer requires a nu- 
merical solution of Brinkman's equations or, equivalently, a solution of the 
(U) problem at high ionic strength. The results in table 1 for uniform layers 
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scaled reciprocal double-layer thickness, Ka scaled reciprocal double-layer thickness, Ka 

Fig. 4. The (scaled) electrophoretic mobility M = 3r]eV/(2e s e kTE) as a function 
of the (scaled) reciprocal double-layer thickness na (aqueous NaCl (I = 1CT 4 -10 M) 
at T = 25°C with radius a « 55,109,219,438,875 and 3500 nm (top-left to 
bottom-right)): a p l 2 = 0.072 with I = 0.71 nm; N = 100; a s = 0.175 A; 
a c ~ —1.95 fiCcm^ 2 . The particles have uniform neutral coatings with thick- 
nesses L and densities uq specified according to Eqn. (60) with F s = 1 and 
Lh ~ (17.5/32)2" nm, where v = 1, 2, 3,..., 8 (labeled in the first panel). 5oZid lines 
are the full model with uniform segment density distributions, and the dash-dotted 
lines are Ohshima's theory (0, Eqn. (11.4.24)). 
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scaled reciprocal double-layer thickness, Ka scaled reciprocal double-layer thickness, Ka 




scaled reciprocal double-layer thickness, Ka scaled reciprocal double-layer thickness, Ka 



Fig. 5. As in figure 4, but with logarithmically scaled mobility axes, 
demonstrate that the hydrodynamic thickness 

L h « L - £, 

where £~ 2 = Qna s F s n . It follows from Eqn. (58) that 

(U - L? ~ ( 1 + L / a ^- 1 
v h ' 18ira s F s Na p /a 
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which is a convenient and useful formula 3 relating the thickness L of a uniform 
layer to the physical parameters L^, Na p , a s and a. The approximation is 
accurate when the 0(j]ME/£) viscous force (per unit volume) on the fluid in 
the layer is small relative to the 0(r]MEL/£ 2 ) frictional (Darcy) drag force, 
i.e., when £/L < 1. 

Equation (60) can be used to reduce the number of unknown degrees of free- 
dom when interpreting electrokinetic experiments. Consider, for example, Co- 
hen and Khorosheva's experiments with poly(ethylene glycol) (PEG) termi- 
nally anchored to the surface of multilamellar liposomes with a m 1750 nm (f2oT). 
Because the PEG grafting density, surface charge density, and PEG molecular 
weight are 'known', Cohen and Khorosheva devised a very simple electroki- 
netic model to infer a value of L and I for each of the five molecular weights of 
PEG, i.e., ten 'fitting' parameters. However, adopting Eqn. (60) reduces the 
number of such parameters to six: one value of L for each molecular weight, 
and one value of a s . 

Recently, Hill interpreted Cohen and Khorosheva's experiments using a self- 



consistent mean-field model of the polymer segment density distribution (Jl 



remarkably small Stokes radius established for PEG statistical segments with 
(Kuhn) length I = 0.71 nm mimics the (negative) correlations that exist be- 
tween micro-scale segment density and fluid velocity fluctuations, which are 
not explicitly accounted for in the Brinkman model (|13l). Physically, regions of 
low segment density favor high fluid velocities, thereby reducing the average 
drag force on an otherwise statistically homogeneous porous medium. This 
suggests that the Stokes radius may be a complicated function of the poly- 
mer structure and density. In relatively dense cross-linked gels, for example, 
density-velocity correlations may be much weaker than in brush-like layers of 
relatively loosely confined polymer. This would lead to a Stokes radius that is 
much more representative of the segment physical size. 

To demonstrate the application of Eqn. (60), we adopt hydrodynamic thick- 
nesses L h and the segment Stokes radius a s ~ 0.175 A as established by 
Hill's self-consistent mean-field description of the polymer (flih. and compare 
the resulting mobilities with Cohen and Khorosheva's experimental data. As 
shown in figure 6, the electrophoretic mobilities predicted in this manner are 
in good agreement with the measured values. Table 3 summarizes the param- 
eters for each polymer molecular weight. Note that all the coatings have the 
same grafting density. Comparing the actual hydrodynamic thicknesses of the 
uniform layers (column 6) with the desired values (column 2) demonstrates 
that the approximation underlying Eqn. (60) is accurate when £/L <C 1, i.e., 
when the PEG molecular weight is greater than 1 kgmol -1 at this grafting 



3 Eqn. (60) is easily solved iteratively, beginning with = L. 




This reduced the problem to one of determining a s 



0.175 A alone. The 
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Table 3 

Parameters characterizing the polymer layers in Cohen and Khorosheva's (CK) 
experiments with terminally anchored PEG on multilamellar liposomes (20). The 
hydrodynamic thickness of the coatings Lh and the Stokes radius of the segments 
a s ~ 0.175 A are specified from Hill's calculations with a self-consistent mean-field 
description of the polymer (flsT). These parameters characterize uniform layers with 
thickness L and permeability £ 2 , specified according to Eqn. (60). 



W (gmol x ) 


N 


L h 1 (nm) 


L (nm) 


n (M) 


I (nm) 


L h 2 (nm) 


350 


4.5 


0.42 


5.5 


0.19 


5.1 


1.5 


1000 


14 


2.0 


4.7 


0.71 


2.7 


2.2 


2000 


28 


4.3 


6.5 


1.0 


2.2 


4.3 


3000 


42 


6.5 


8.6 


1.2 


2.1 


6.5 


5000 


70 


11 


13 


1.3 


2.0 


11 



Actual or desired values from Hill's self-consistent mean-field calculations (II 



2 Actual values for uniform layers. 

density. The results are also consistent with expectations that the density 
and, hence, permeability of polymer brushes are independent of the molecu- 
lar weight at a fixed grafting density. Clearly, the Brinkman screening length 
decreases with increasing molecular weight, approaching I « 2 nm. The layer 
thickness increases linearly with molecular weight, yielding L « (0.21+0. 22N)l 
(I ~ 0.71 nm) when W > 1 kgmol -1 . As expected, this formula is the same as 
Hill's for the high-molecular-weight limit of the hydrodynamic layer thickness, 
Lh ~ 0.22NI, at this grafting density. 

3.5 Terminally anchored PEG on lipid bilayer membranes (stealth liposomes) 

Cohen and Khorosheva's experiments, introduced above, were performed with 
relatively large multilamellar liposomes. This permitted the particles to be 
observed under magnification in their electrophoresis device. However, stealth 
liposomes used for drug delivery have much smaller radii a ~ 100 nm. With 
the same surface charge, polymer, and polymer grafting density, the smaller 
radius is expected to change the surface potential and, possibly, vary the 
polymer segment density distribution. Therefore, with knowledge of the Stokes 
radius for PEG (Kuhn) segments, obtained from Cohen and Khorosheva's 
experiments with a = 1750 nm, we can predict the electrophoretic mobilities 
of therapeutic stealth liposomes with a = 100 nm, for example. 

The polymer segment density distributions (figure 7) are from Hill's self- 
consistent mean- field model (|13T). These calculations account for surface cur- 
vature, yielding slightly less dense layers with a = 100 nm than with a = 
1750 nm. The electrostatic potential is much lower, and, as expected, this 
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scaled recipocal double-layer thickness, km 



Fig. 6. The (scaled) electrophoretic mobility M = 3neV/(2e s e kTE) of spher- 
ical liposomes with coatings of terminally anchored PEG as a function of the 
(scaled) reciprocal double-layer thickness na (aqueous NaCl at T = 25° C with 
radius a = 1750 nm) for various numbers of statistical segments per PEG chain 
N = 0, 4.5, 14, 28, 42, and 70 (increasing downward); the molecular weights of 
the PEG chains are, respectively, W = 0, 0.35, 1, 2, 3 and 5 kgmol -1 . The surface 
charge density a c ~ —1.95 /xCcm -2 at all ionic strengths. Circles (with dotted lines 
to guide the eye) are Cohen and Khorosheva's experimental data ([20]) . and solid 
lines are the full model with uniform segment density distributions, whose density 
and thickness are chosen to yield the same hydrodynamic layer thicknesses as Hill's 
self-consistent mean-field description of the layers |l3h. See table 3 for parameters. 

changes the electrophoretic mobility, as seen by comparing the mobilities in 
figures 6 and 8 with a = 1750 and 100 nm, respectively. 



3. 6 PEO adsorbed on polystyrene latices 



As demonstrated in section 3.4, the brush-like structure of terminally anchored 
PEG permits a satisfactory uniform-layer approximation when the density and 
thickness yield the 'correct' hydrodynamic coating thickness. Poly (ethylene 
oxide) (PEO) homopolymer, which comprises the same monomer unit as PEG, 
has molecular weights an order of magnitude higher. Adsorption produces 
much more inhomogeneous layers (see 21 ; 2^, and the references therein) , and 



because desorption is slow, experiments can be performed under conditions 
where the influence of polymer on the hydrodynamic size and electrophoretic 
mobility is significant, without altering the electrolyte viscosity. 



Gittings and Saville (23) measured the adsorbed amounts, hydrodynamic layer 
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scaled distance, (r-a)/\ 



Fig. 7. The (scaled) radial segment density distribution (ft = n(r)/ 3 of PEG chains 
terminally anchored to the surface of spherical (therapeutic) liposomes with ra- 
dius a = 100 nm as a function of the (scaled) distance from the surface of the 
bare liposome (r — a) /I. The data are for representative numbers of statistical seg- 
ments per chain N = 28, 70, 140 and 280, with molecular weights W 2, 5, 10 
and 20 kgmol -1 . The self-consistent mean-field potential is specified according to 
Hill (Q) with wa/ui = 3.3, v 2 /l 3 = 0.27, a p l 2 = 0.072 and I = 0.71 nm. 

thicknesses, and electrophoretic mobilities of polystyrene latices in aqueous 
solution with three molecular weights of adsorbed PEO. They used the mo- 
bilities of the bare particles to infer the (varying) surface charge density as 
a function of ionic strength. As shown in table 4, the (apparent) charge, in- 
ferred from the standard electrokinetic model and measured mobilities of bare 
latices, increases with ionic strength. 

With the methodology presented in section 3.4, these very inhomogeneous 
structures are approximated as uniform layers. Also, following Gittings and 
Saville, this analysis assumes that the polymer does not influence the surface 
charge. Table 5 summarizes the measured and inferred layer characteristics, 
and the electrophoretic mobilities and effective coating thicknesses are shown 
in figures 9 and 10. 

Because adsorbed PEO produces layers with a very dense inner region and 
a much more permeable periphery, approximating the layers as uniform is 
crude. Nevertheless, varying the segment density distributions to improve the 
predicted mobilities is invariably met with a poorer prediction of the measured 
effective layer thicknesses. In particular, adopting a model with two step-like 
regions (a very dense and relatively thin inner layer with a much thicker and 
permeable outer layer) produced results very similar to those in figures 9(a) 
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scaled reciprocal double-layer thickness, Ka 



Fig. 8. The (scaled) electrophoretic mobility M = 3rjeV/(2e s e kTE) of spher- 
ical liposomes with coatings of terminally anchored PEG as a function of the 
(scaled) reciprocal double-layer thickness Ka (aqueous NaCl at T = 25° C with 
radius a = 100 nm) for various numbers of statistical segments per PEG chain 
N = 4.5, 14, 28, 70, 140 and 280; the molecular weights of the PEG chains are, 
respectively, W « 0.35, 1, 2, 5, 10 and 20 kgmol -1 . The surface charge density 
a c ~ —1.95 [iCcmT 2 at all ionic strengths. These predictions are from the full 
model with a self-consistent mean-field description of the PEG segment density 
distribution (see figure 7). 

and 10(a) with a single region. This was also the case with exponentially de- 
caying distributions. Varying the Stokes radius of the segments, while main- 
taining the 'correct' hydrodynamic layer thicknesses and adsorbed amounts, 
slightly improved the predicted electrophoretic mobilities, but the resulting 
layer thicknesses L were unrealistically large, and the resulting Stokes radius 
a s was extremely small. 

Comparing the model predictions in figure 9(a) with Gittings and Saville's 
measurements in figure 9(b) reveals that the influence of polymer on the elec- 
trophoretic mobilities is underestimated by the model at low ionic strength, 
and overestimated at high ionic strength. The qualitative effects of varying 
the PEO molecular weight and adsorbed amount are captured relatively well, 
however. Comparing variations in the effective layer thicknesses of the bare 
and coated particles in figure 10 demonstrates that polymer attenuates the 
electro-viscous contribution to the drag. Even for the bare latices, variations 
in the effective size with ionic strength are greater than suggested by the the- 
ory. Recall, L e reflects an increase in the apparent radius of the bare particle. 
Therefore, for colloids without a polymer coating, L e equals the hydrodynamic 
thickness Lh only in the absence of electro-viscous drag. 
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Table 4 

The ionic strength-charge relation for the PEO-coated polystyrene latices (a = 
78 nm) used in the Gittings-Saville experiments (23). The surface charge densities 
are inferred from the measured electrophoretic mobilities of the bare latices. See 
table 5 for the PEO-coating parameters. 

/ (M) a c (/iCcm~ 2 ) 

1(T 5 0.250 

10~ 4 0.600 

5 x 10" 4 0.995 

10~ 3 1.352 

3.43 x 10~ 3 1.713 

10~ 2 1.928 



Evidently, the influence of adsorbed polymer on the surface charge depends 
on the ionic strength. In general, therefore, neutral polymer tends to reduce 
the (apparent) charge more than expected from the effect on hydrodynamic 
transport processes alone. Either the polymer changes the underlying charge 
or charge affects the polymer configuration. At very low ionic strength, a con- 
siderable portion of the diffuse double layer resides outside a uniform polymer 
layer, suggesting that either the polymer extends further from the interface, 
to increase the effective layer thickness, or the charge is diminished. 



Note that polymer lowers the effective dielectric constant at the interface, 
thereby increasing correlations between the fixed and mobile charge, which, 
in turn, could decrease the degree of counterion dissociation. Note also that 
polarization of low-dielectric neutral polymer segments by the equilibrium 
electrostatic potential may increase layer thicknesses, particularly at low ionic 
strength when the double layer is thick and the surface potential is high. 
We have not explored the possibility that the Stokes radius of the segments 
may vary with the local polymer configuration. With large fluctuations in seg- 
ment density in the periphery, the effective Stokes radius might be very small, 
whereas the thin dense region at the surface may produce an effective Stokes 
radius more representative of the segment size. The model also neglects the 
influence that dense polymer has on molecular diffusion and electro-migration. 
Therefore, while the full model improves upon earlier interpretations based on 
the standard model (see Eih . further investigation is required. 
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(a) (b) 

Fig. 9. Theoretical interpretation (left) of Gittings and Saville's experiments (right) 
for the electrophoretic mobility V/E of polystyrene latices with coatings of adsorbed 
PEO as a function of the ionic strength / (aqueous NaCl at T = 25° C) for various 
PEO molecular weights W = (open circles), 23.5 (filled circles), 56 (squares) and 
93.75 kgmol -1 (triangles). The surface charge densities and adsorbed amounts are 
specified according to Gittings and Saville (J23J) (see tables 5 and 4), and the mobil- 
ities are calculated using the full electrokinetic model with uniform layers yielding 
the measured hydrodynamic layer thicknesses and adsorbed amounts. The segment 
Stokes radius a s = 0.175 A with I = 0.71 nm is from Hill's interpretation |l3h 
of Cohen and Khorosheva's experiments (|2ol ) with terminally anchored PEG on 
multilamellar liposomes. 



Table 5 

Polymer layer properties for the PEO-coated polystyrene latices used in the 
Gittings- Saville experiments |23j). The electrophoretic mobilities of the bare and 
'fuzzy' latices are presented in figure 9. These parameters characterize the uniform 
layers adopted in the full electrokinetic model. See table 4 for the surface charge 
densities. 



W (kgmol x ) 


cr p (mgm -2 ) 


L h 1 (nm) 


N 


L (nm) 


£ (nm) 


L h 2 (nm) 


23 


0.67 


11.5 


534 


13.8 


2.32 


11.4 


56 


0.73 


17.0 


1270 


19.8 


2.76 


16.9 


93.75 


0.59 


23.5 


2130 


27.3 


3.76 


23.3 



1 From Gittings and Saville's experiments (j23j). 

2 Calculated with uniform layers. 
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Fig. 10. Theoretical interpretation (left) of Gittings and Saville's experiments (right) 
for the effective layer thickness L e of PEO adsorbed on polystyrene latices as a func- 
tion of the ionic strength I (aqueous NaCl at T = 25° C) for various PEO molecular 
weights W = (open circles), 23.5 (filled circles), 56 (squares) and 93.75 kgmol -1 
(triangles). The surface charge densities and adsorbed amounts are specified accord- 
ing to Gittings and Saville (|23T ) (see tables 5 and 4), and the effective layer thick- 
nesses are calculated using the full electrokinetic model with uniform layers yield- 
ing the (nominal) measured hydrodynamic layer thicknesses and adsorbed amounts. 
The se gme nt Stokes radius a s = 0.175 A with I = 0.71 nm is from Hill's interpre- 
tation (|l3h of Cohen and Khorosheva's experiments |2oh with terminally anchored 
PEG on multilamellar liposomes. Lines indicate the hydrodynamic layer thicknesses 



4 The mobility of colloids with charged coatings 

4-1 Thin poly electrolyte coatings 

In this section we examine moderately large particles with thin charged (poly- 
electrolyte) layers. The electrophoretic mobilities of colloids with radius a = 
500 nm in figure 11 show the effect of varying the (uniform) polyelectrolyte 
density while maintaining a constant layer thickness L = 0.25-\/l0 3 nm. As 
the ionic strength increases from / = 10~ 6 to 10 M, na spans the range 1.64- 
5200; the double-layer thickness is comparable to the coating thickness when 
/ pa 1.5 mM. The densities of the coatings span the range Uq pa 0.13-1.1 M, 
so with a segment Stokes radius a s = 0.95 A the respective Brinkman screen- 
ing lengths span the range £ pa 2.6-0.89 nm. Note that Stokes radius of the 
(monomer) segments adopted here is an approximate value for poly(styrene 
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sulfonate) monomer segments inferred by Hill (|13T ) from experiments by Cot- 
tet et al. (jl^) reporting the electrophoretic mobilities and hydrodynamic radii 
of associating copolymer micelles (|13l ) (see § 4.3). 

The fixed charge density is assumed proportional to the segment (monomer) 
density and, furthermore, is specified assuming one elementary charge per 
Bjerrum length Ib of the polyelectrolyte contour length. Accordingly, the 
charge density nl = (2.5/7.1)no; each segment has length I = 2.5 A with 
Ib = e 2 / (4iie e s kT) pa 7.1 A. This simple model accounts for counterion con- 
densation, which, according to Manning's well-known theory (j25h. limits the 
effective linear charge density to e/^ 1 when the actual linear charge density 
exceeds this value. For highly charged polyelectrolytes, e.g., poly(styrene sul- 
fonate) with degree of sulfonation greater than 40 percent, accounting for 
counterion condensation in this way leads to 'reasonable' electrostatic poten- 
tials (Donnan potentials) inside the layer. When interpreting measured elec- 
trophoretic mobilities, the methodology infers a Stokes radius for the monomer 
segments a s pa 0.95 A, which, in contrast to PEG brushes (see § 3.4 and § 3.5), 
is representative of the segment (monomer) length, l m pa 2.5 A. 

The moderately large particle radius (a = 500 nm) facilitates a comparison 
of the full model with analytical solutions valid as na — > oo with kL ^> 1. 
Levine et al.'s theory (0, Eqn. (21)) (dashed lines) is in excellent agreement 
with the numerically exact results when the ionic strength is high, i.e., when 
kL ^> 1 and the Debye-Hiickel approximation is justified (ip < 2kT/e). The 
range of ionic strengths over which their formula is accurate increases with 
decreasing charge density, suggesting that the Debye-Hiickel approximation 
is more restrictive here than the condition kL ^> 1. Ohshima's theory (flil 
Eqn. (11.4.18)) (dash-dotted lines), which neglects the viscous stress at the 
grafting surface, is restricted here to dense layers with L/£ > 8. 

At high ionic strength, viscous drag at the grafting surface leads to lower mo- 
bilities than predicted by Ohshima's theory. This is demonstrated most clearly 
by the particle with the least dense layer 4 (case 1 in figure 11). Under free- 
draining conditions — the regime that underlies the high-ionic-strength-limit 
of Ohshima's theory — the electrophoretic mobility of this particle is expected 
to be the highest, since the drag coefficient of the segments is least because of 
the small, but finite, (hydrodynamic) volume fraction <p s = n(4/3)nal ~ 10~ 4 . 
However, because the layer is sufficiently permeable and thin, viscous stress 
from the substrate also contributes to the drag. Consequently, this particle 
has the lowest mobility at high ionic strength. At lower (intermediate) ionic 
strength, an increasing fraction of the mobile charge resides beyond the coat- 



Close examination of figure 11 reveals that Ohshima's theory (dash-dotted lines) 
yields a mobility about 10 percent higher (in magnitude) than the full theory at 
high ionic strength. 
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Fig. 11. The (scaled) electrophoretic mobility M = 3rjeV/(2e s e kTE) of spherical 
colloids with relatively thin, uniform coatings of charged polymer as a function of 
the ionic strength / (aqueous NaCl at T = 25°C): o = 500 nm; L = 0.25V 10 3 nm; 
a s = 0.95 A; n « 0.132, 0.263, 0.526 and 1.05 M (labeled 1-4 with i « 2.6, 1.8, 
1.3 and 0.89 nm); and = (2.5/7.1)no- Solid lines are the full model with uniform 
(step-like) segment density distributions, dash- dotted lines are Ohshima's theory 
for thin uniform layers with L/i ^> 1 (J19J, Eqn. (11.4.18)), and dashed lines are 
Levine et al.'s theory for thin uniform layers with arbitrary L/i and < 2kT/e (0, 
Eqn. (21)). 

ing, and the mobility might therefore be expected to increase with the total 
charge (increasing segment density). However, the particles with the highest 
charge are more susceptible to polarization, and this leads to an inverse re- 
lationship between the charge and mobility. Nevertheless, at even lower ionic 
strength, particles with higher charge eventually exhibit higher mobilities. In 
this limit, the total charge, rather than the slightly larger hydrodynamic ra- 
dius (friction coefficient), has the greater influence on the mobility. 



4-2 Glycocaylx on human erythrocytes 

Levine et al. (5) developed their theory to interpret the electrophoretic mo- 
bility of human erythrocytes, revealing that the glycocaylx charge density is 
sufficiently high to produce an electrostaticpotential that violates the Debye- 
Hiickel approximation. Sharp and Brooks (p) then solved the planar electroki- 
netic transport problem numerically, relaxing the Debye-Huckel approxima- 
tion, and they were able to infer a layer thickness of L = 7.8 A, with segment 
and charge densities of n ~ 0.0690 and nl w 0.0472 M, respectively, at ionic 
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strengths in the range I = 10-160 mM. Note that, with a segment Stokes 
radius a s = 7 A and drag coefficient F s — 1, the Brinkman screening length is 
£ w 1.35 nm. Sharp and Brooks' calculations are limited by the finite particle 
size (polarization and relaxation), so it remains to establish the range of ionic 
strengths over which the flat-plate approximation is accurate. We will also 
briefly examine the effect of varying the layer thickness, while maintaining a 
constant mass and charge. 

The mobility-ionic strength relationship in figure 12 verifies the accuracy of 
Ohshima's theory (dash-dotted lines) for this problem when the ionic strength 
is greater than approximately 5 mM. Clearly, Levine et al.'s theory (dashed 
lines) is limited here by the Debye-Hiickel approximation, even at low ionic 
strength. At high ionic strength, however, a close examination reveals that 
their expression is more accurate than Ohshima's, because it accounts for 
viscous stress on the bare surface. Finally, the 'exact' results (solid lines) show 
that the particle size (a = 3.5 /jm) influences the mobility when the ionic 
strength is less than about 10 mM, i.e., when na < 10 3 . 

To assess the effect of expansion and contraction of the layers due to elec- 
trostatic stiffening of the polyelectrolyte at low ionic strength, we perturbed 
Sharp and Brooks' data (case 2) by halving (case 1) and doubling (case 3) the 
(glycocaylx) layer thickness. Because the particle radius is large (a = 3.5 /im), 
the changes are accompanied by respective doubling and halving of the seg- 
ment and charge densities. Expansion decreases the segment density, thereby 
decreasing the fixed charge density and increasing the layer permeability. To 
account for the influence of density on the segment drag coefficient, the Stokes 
radius of the segments is a s = 3.5 A, which yields a slightly larger Brinkman 
screening length (£ = 1.73 A) than that of Sharp and Brooks when L = 7.8 A. 
The hydrodynamic volume fraction inferred by Sharp and Brooks' parameters 
suggests that varying the density should also affect the segment drag coeffi- 
cient. Indeed, we have interpreted their Stokes radius as the product a s F s (<j) s ), 
so changing the density also changes the drag coefficient according to the 
change in (hydrodynamic) volume fraction. The sensitivity of the mobility to 
the layer thickness and, hence, density, is most pronounced at an ionic strength 
of approximately 1 mM. Because polarization and relaxation are significant 
under these conditions, numerically exact solutions could benefit future stud- 
ies, even with relatively large particles. 



4-3 Poly(styrene sulfonate) micelles 

Finally we turn to 'small' colloids with relatively thick charged (sodium poly(styrene 
sulfonate)) layers. Figure 13 shows the mobilities of two representative hy- 
drophobically associating copolymer micelles over a range of ionic strengths. 
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Fig. 12. The electrophoretic mobility V/E of human erythrocytes as a function 
of the ionic strength / (aqueous NaCl at T = 25°C): a = 3.5 ^m; L = 7.8/2, 7.8 
and 7.8 x 2 nm (labeled 1-3, respectively); a s = 3.5 A; n ^ 0.0690 x 2, 0.0690 
and 0.0690/2 M (labeled 1-3 with £ 2.52, 1.73 and 1.16 nm, respectively); and 
71q = (0. 0472/0. 0690)no . Solid lines are the full model with uniform (step-like) seg- 
ment density distributions, dash-dotted lines are Ohshima's theory for thin uniform 
layers with L/i ^> 1 (|13 . Eqn. (11.4.18)), and dashed lines are Levine et al.'s theory 
for thin uniform layers with arbitrary L/£ and < 2kT/e 0, Eqn. (21)). 

The radii of the hydrophobic cores, a, and of the respective particles as a whole, 
Lh + a, are based on neutron diffusion and light scattering measurements, re- 
spectively (|24f ). Furthermore, the charged coronas are assumed uniform, with 
monomer densities 



n 



n N 
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(61) 



where n a is the aggregation number and N m is the number of poly(styrene 
sulfonate) monomers per chain. Note that L and Lh are specified according to 
Eqn. (60); these and other characteristics are listed in table 6. 

The micelles assembled from the longest polyeiectrolyte chains (case 2) have 
a much larger hydrodynamic radius + a w 80 nm. The smaller aggregation 
number n a pd 35 manifests in a relatively low segment density n Q w 0.0074 M, 
and, hence, lower density of charged sites nl w (2.5/7. l)n w 0.0026 M. Note 
that, following the discussion in section 4.1, allowing for counter ion conden- 
sation yields an effective valence for monomer segments z ~ —2.5/7.1. The 
total charge, which is proportional to zN m n a , is only 30 percent lower, how- 
ever. Clearly, the lower segment density results in a much more permeable 
particle, as evidenced by the large Brinkman screening length i w 11.2 nm. 
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Fig. 13. The (scaled) electrophoretic mobility M = 3r]eV/(2e s e kTE) of spherical 
colloids with relatively thick, uniform coatings of charged polymer as a function of 
the ionic strength / (aqueous NaCl at T = 25° C). See table 6 for parameters. For 
both particles, a s = 0.95 A and 4 = (2.5/7.1)no- Solid lines are the full model 
with uniform (step-like) segment density distributions, and dash-dotted lines are 
Ohshima's theory for spherical polyelectrolytes Eqn. (11.3.11)). 

Table 6 

Parameters derived from data (|24h characterizing two micelles (with charged coro- 
nas) whose electrophoretic mobilities are shown in figure 13: a s = 0.95 A and 
4 = (2.5/7.1)7*,. 



case 


a (nm) 


N 


n a 


Nn a 10- 3 


Lh + a 1 (nm) 


1 


5.2 


227 


85 


19.3 


42 


2 


2.6 


404 


35 


14.1 


80 



case 


L + a 2 (nm) 


n (M) 


4 ( M ) 


£ (nm) 


L/£ 


L/a 


1 


45.4 


0.082 


0.029 


3.31 


12 


7.7 


2 


91.2 


0.0074 


0.0026 


11.2 


7.9 
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1 Actual (measured) hydrodynamic size. 

2 From the measured hydrodynamic size and Eqn. (60). 



Following the methodology underlying Eqn. (60), the actual 'coating' thick- 
ness L 89 nm is about 12 percent smaller than the chain contour length 
l c m 404 x 0.25 = 101 nm, indicating that the chains are strongly stretched. 

Note that the permeabilities in table 6 have been established from electrophoretic 
mobilities reported by Cottet et al. (24) (at a single ionic strength yielding 



32 



k- 1 ta 2 nm): V/E « -3.4 x 10~ 8 (case 1) and -3.7 x 10~ 8 m 2 /(Vs) (case 2). 
From the high-ionic-strength limit of the full electrokinetic model 

V/E = ze/{67rr]a s F s ), (62) 

where z is the effective valence of a segment, these mobilities and the 'coating' 
parameters in table 6 suggest a Stokes radius for monomer segments a s « 
0.95 A (Q. 

Equation (62) is the free-draining limit of Ohshima's (flil Eqn. (11.3.11)) 
and Hermans and Fujita's theories for uniformly charged spherical polyelec- 
trolytes (19). Solutions of the full electrokinetic model show that Eqn. (62) 
also applies to non-uniform layers when the segment and charge densities are 
proportional to each other, and the coatings are thick enough for the viscous 
stress on the substrate to be small compared to the frictional drag on the 
polymer fl^K 

Figure 13 shows that the mobilities predicted by the full model do not vary 
significantly over a wide range of electrolyte concentrations. While Ohshima's 
theory (dash- dotted lines) is in good agreement with the 'exact' results at 
very high ionic strength, the correction it provides to Eqn. (62) is very small. 
The analytical theory compares more favorably with the 'exact' results for 
the larger particle, i.e., case 1 with larger n(a + L), despite the smaller value 
of L/i. Again, this points to the significance of polarization and relaxation. 
Note that the high ionic strength at which Eqn. (62) applies usually permits 
polarization and relaxation to be neglected. For these relatively small and 
highly charged particles, however, the analytical theory is accurate only at 
inordinately high ionic strength. 

In contrast to bare particles and, indeed, particles with neutral coatings, the 
mobilities reach minimum values at intermediate ionic strengths. Here, the 
minima coincide with the electrostatic potential inside the polyelectrolyte 
passing through Pd 2kT/e. Note that the potentials would be significantly 
higher if we had not accounted for counterion condensation, at least in the 
manner described in section 4.1. With a characteristic particle radius a + L^, 
n(Lh + a) ~ 14 (case 1) and 26 (case 2) when / = 10 mM, for example. More- 
over, when / = 1 M, n(Lh + a) « 140 (case 1) and 260 (case 2). Once again, 
with the mobilities shown in figure 13, we see that polarization and relaxation 
is important when the double-layer thickness is greater than a few percent of 
the characteristic particle size. 

At low ionic strength, the mobility reflects a balance between the electrical 
force on the fixed charge, which is proportional to zNn a , and the hydrody- 
namic drag, which is proportional to a + Lh ~ N u (0 < v < 1). It follows 
that the mobility should increase with the aggregation number, which is in- 
deed the case at low ionic strength. Because the mobility is also susceptible to 
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polarization, and not all the counter charge resides beyond the polymer, the 
permeability and particle size also play a role. 

As expected, calculations (not reported here) with twice the segment Stokes 
radius (a s = 1.9 A) yield mobilities (at high ionic strength) that are half 
those in figure 13. At low ionic strength, the mobilities are attenuated only 
slightly, whereas at intermediate ionic strengths the minima cease to exist 
and the mobility (magnitude) increases monotonically with decreasing ionic 
strength. Evidently, decreasing the permeability decreases the effectiveness 
of convection in polarizing the double layer, as expected from the inferences 
drawn in section 3.3 for neutral layers. 

Note that the thickness and density of terminally anchored polyelectrolytes 
can vary significantly with the bulk ionic strength. Such variations have been 
neglected entirely. At low ionic strength, the effective Kuhn length increases 
because of weakly screened electrostatic repulsion among charged monomers 
along the backbone (a so-called short-range interaction). This can lead to 
layers whose hydrodynamic thickness is comparable to the polymer contour 
length. Here, however, the relatively high ionic strength at which the mobilities 
and sizes were reported suggests that excluded volume (long-range) interac- 
tions are dominant because of the high grafting density. At lower densities, 
electrostatic interactions can be much more influential (see Eil ). Moreover, 
grafted chains may 'collapse' at high ionic strength if the grafting density 
is low, yielding much thinner, dense layers. The varying segment and charge 
density that accompany such structural changes also affect the electrostatic 
potential, which, in turn, influences the ionic strength at which polarization 
affects the mobility. A quantitative interpretation of these processes would 
benefit greatly from experiments reporting mobilities and hydrodynamic radii 
over a range of ionic strengths. 



5 Summary 

The 'full' electrokinetic model of Hill, Saville and Russel (flih was applied to a 
variety of 'soft' colloids. Analogously to O'Brien and White's solutions of the 
standard electrokinetic model, the methodology removes all approximations 
imposed by earlier theories. Because of polarization and relaxation, exact solu- 
tions of the full model were demonstrated to be as important for 'soft' colloids 
as the O'Brien and White methodology is for 'bare' particles. 

A simple approach was described to link characteristics, such as the poly- 
mer adsorbed amount and molecular weight, to key parameters in the model, 
namely the polymer layer density, permeability and thickness, and their radial 
distributions. Calculations revealed that the hydrodynamic layer thickness, as 
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distinguished from the actual thickness, is the single most important layer 
characteristic influencing the mobility. A convenient (and often reasonable) 
approximation is to set the hydrodynamic thickness of a uniform layer to its 
actual thickness minus the Brinkman screening length. It is then straightfor- 
ward (e.g., Eqn. 60) to relate the actual thickness, density and permeability 
(via the Stokes radius of the segments) of a uniform coating to the grafting 
density (or adsorbed amount) and polymer molecular weight for a non-uniform 
layer. 

The full model was compared to Ohshima's analytical theories for thin uni- 
form, neutral layers. As expected, good agreement was found when the double 
layer and coating are very thin compared to the particle radius (less than a 
few percent) and polarization and relaxation are negligible. Figures 4 and 5 
provide a useful reference for assessing the (coupled) influences of layer thick- 
ness, double-layer thickness and particle size on polarization and relaxation. 
Contrary to expectations, flat-plate approximations breakdown at very high 
ionic strength. However, because the mobility under these conditions is small, 
this unexpected behavior is revealed only on logarithmically scaled mobility 
axes. Neutral polymer was identified as attenuating polarization, but polar- 
ization and relaxation still influence the mobility of particles with radii less 
than a few microns. 

The full model was applied to interpret the mobilities of (stealth) liposomes 
with terminally anchored poly (ethylene glycol) (PEG). Approximating the 
layers as uniform, with the 'correct' hydrodynamic layer thickness, provided 
a close correspondence with experiments (J20J) and calculations that invoke 
much more involved computations of the polymer- layer structure ffl3h . There 
remain relatively small differences between theory and experiment when the 
ionic strength and polymer grafting density are low ffl3h . Polymer anchored to 
the surface of fluid-like membranes poses complications beyond the scope of 
the current elect rokinetic model, which assumes immobile surface charge and 
polymer. Nevertheless, based on earlier attempts to establish the Stokes radius 
of aqueous PEG statistical segments, theoretical predictions of the mobilities 
of 100 nm radius (therapeutic) stealth liposomes were presented with a typical 
grafting density and a variety of polymer molecular weights. These calculations 
clearly demonstrate the extent to which neutral polymer coatings decrease the 
apparent ^-potential, even though the actual charge remains constant. 

When applied to interpret the mobilities and hydrodynamic sizes of polystyrene 
latices stabilized in aqueous electrolyte by adsorbed poly(ethylene oxide) (PEO) 
homopolymer, the theory captured only the qualitative aspects of the avail- 
able data Q. Moreover, measured mobilities at low ionic strength are lower 
than suggested by the theory. Evidently, adsorption of neutral polymer lowers 
the charge, or the electric field due to charge alters the polymer conformation. 
While model parameters (for the polymer layers) can be adjusted to provide a 
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good fit to measured electrophoretic mobilities, this is inevitably accompanied 
by a poorer prediction of the (measured) hydrodynamic layer thickness. 

For particles with charged layers, the full model corroborated approximate 
theories in their respective domains of applicability (e.g., thin layers, at high 
ionic strength, with low electrostatic potentials). For example, we adopted 
parameters that characterize human erythrocytes, as established from theo- 
retical and experimental studies by Brooks and coworkers. We also compared 
the full model with approximate analytical solutions for particles with 'thick' 
charged layers. These parameters were based on measured dimensions (at one 
ionic strength) of polyelectrolyte micelles (with poly(styrene sulfonate) coro- 
nas). Extrapolating the electrophoretic mobilities to higher and lower ionic 
strength (with fixed dimensions and charge) revealed that analytical theory 
(after Ohshima and Hermans and Fujita) is accurate (for these particles) only 
at inordinately high ionic strength. This shortcoming reflects the small size 
of the micelles, their high charge density, and, hence, the strong influence of 
polarization and relaxation. Note that the calculations invoked counterion con- 
densation, which places an upper limit on the effective charge and, therefore, 
limits the prevailing electrostatic potential. 

At present, the availability of sufficiently well-characterized electrokinetic data 
for soft spheres is limited, so we trust that the model and its availability will 
stimulate further experimental studies. Again, it should be noted that, while 
this paper focused on the electrophoretic mobility, the electrokinetic model 
also provides characteristics (e.g., complex polarizability/conductivity (jl^lH) 
and dynamic mobility) derived from a variety of electrokinetic experiments. 

RJH acknowledges support from the Natural Sciences and Engineering Re- 
search Council of Canada (NSERC), through grant number 204542, and the 
Canada Research Chairs program (Tier II). 
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